Quantum System Identification 
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The aim of quantum system identification is to estimate the ingredients inside a black box, in 
which some quantum-mechanical unitary process takes place, by just looking at its input-output 
behavior. Here we establish a basic and general framework for quantum system identification, that 
allows us to classify how much knowledge about the quantum system is attainable, in principle, from 
a given experimental setup. Prior knowledge on some elements of the black box helps the system 
identification. We present an example in which a Bell measurement is more efficient to identify the 
system. When the topology of the system is known, the framework enables us to establish a general 
criterion for the estimability of the coupling constants in its Hamiltonian. 



Some of the most exciting and puzzling concepts in 
quantum theory can already be observed in simple sys- 
tems. These are, for example, superpositions and deco- 
herence, tunneling, entanglement and non-locality, quan- 
tum cryptography, teleportation and dense coding. Many 
of such theoretical ideas have been confirmed experimen- 
tally with a tremendous accuracy. On the other hand, 
perhaps the most important theoretical concept — a full 
quantum computer or simulator — is still well out of reach, 
because it requires a fully controllable system of Hilbert 
space dimension at the very least of the order of 2^"°. Its 
realization poses one of the greatest challenges in science 
today. 

On our path towards quantum computation we are 
building systems composed of more and more qubits, 
the quantum information theoretic equivalent of the bit. 
But while an information theoretic approach is very suc- 
cessful, we should not forget that any implementation 
comes with a baggage of physical effects. In particu- 
lar, real qubits interact. Often, these interactions are 
important: they are actively used to create logical gates. 
Sometimes, they are unwanted, and either suppressed ac- 
tively, or simply neglected. However if we are to meet the 
stringent bounds that fault-tolerance computation puts 
on the required precision of our technology, we will have 
to estimate our quantum system with very high preci- 
sion. Current estimates of the fault-tolerance threshold 
indicate that in many systems the relative precision will 
have to be of the order of 10~^-10~^. 

If we could perfectly control our system, achieving such 
precisions is a mere engineering difficulty. But if our con- 
trol relies on the system couplings, or is heavily perturbed 
by them, we are in a Catch-22 situation, and it is unclear 
how well the system can be estimated even in principle. 
In this paper, we solve this question by providing a pre- 
cise mathematical description of the equivalent set [l|-|3| 
of closed systems. This set describes the possible im- 
plementations of a system that cannot be distinguished 
with a given experimental setup. It should be compared 
to the well-known reachable set in quantum control 4], 
which describes the set of unitary operations that can be 
implemented, in principle, by a given experimental setup. 



It has been shown in quantum control that even when 
only parts of the system are accessed, the reachable set 
typically remains maximal: the system is capable of 
quantum computation We show that this is not true 
for full estimability: in general, infinitely many different 
system Hamiltonians give rise to the same input-output 
behavior. However, we show how a priori knowledge 
about the system helps to restrict the set of possible sys- 
tems. Indeed we prove that in a generic limited- access 
situation, relatively little a priori knowledge can imply 
full estimability. This generalizes several recently devel- 
oped schemes for indirect estimation [gI-Q . We also show 
how estimability can strongly depend on the structure of 
quantum measurements, by providing an example where 
entangled observables are more efficient for the estima- 
tion than product observables. 

Our analysis first follows closely the known results from 
bilinear theory [ij. Then, we use a result from Lie alge- 
bras to translate the bilinear theory to the quantum 
case. This sets our result apart from previous work which 
required additional mathematical assumptions 0, Q . 

Setup. — We consider a black box with Ni inputs 
and No outputs. Inside the black box, some quantum- 
mechanical unitary dynamics takes place. Our goal is to 
find a model for the black box that perfectly describes its 
input-output behavior under all possible circumstances 
[system identification [Tl). 
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Figure 1: A set of time dependent functions fk{t) is the 
input, which determines the unitary dynamics inside the 
black box, and a set of the expectation values of observables 
Mi is the output. Our objective is to estimate the system 
a = {Ho, Hk, Me, Po} by looking at the input-output behav- 
ior of the black box. In the most extreme case, even the 
control operations Hk and the observables Mi are unknown. 
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More specifically, we are modeling a system with a 
finite-dimensional Hilbert space H, a time dependent 
Hamiltonian 

Ni 

fc=i 

an initial quantum state po, and a set of observables Mi 
(£ = !,..., No). Without loss of generality we chose Hq 
and Hk traceless. The inputs are the functions /fe(i) (fc = 
1, . . . ,Ni), which are assumed to be piecewise constant. 
The outputs are the expectation values of the observables 
Mi, 

tr{Mip{t)} with p{t) = T^exp(^J dt' C{t')^ po, 
where 

Ni 

C{t)= C^+Y,fk{t)Ck, Ck = -i[Hk,*] (fc = 0,...,iV,) 

k=l 

are the Liouvillians corresponding to the Hamiltonians. 
See Fig. [1] Because we are interested in whether systems 
can be distinguished in principle, we assume that it is 
possible to collect statistics at arbitrary precision, and 
that infinitely many copies of the system are available 
(this allows us to ignore any back-action of the measure- 
ments 0, 0] ) ■ Our main assumption is that the system 
is controllable, implying that any unitary transformation 
can be realized by the Hamiltonian dynamics with H{t), 
by properly arranging the inputs /fc(i). Mathematically 
this amounts to the smallest Lie algebra over the re- 
als that contains the matrices iHq, iHi, . . . , 01^^ being 
equal to the full Lie algebra su(dim?^) of traceless skew- 
hermitian matrices of size dimH x dimH. Controllability 
is a generic property of systems, and is in principle an 
observable property, if the dimension of the underlying 
Hilbert space is known. Furthermore we exclude the triv- 
ial cases where or pQ is proportional to the identity 
operator. 

We put all parameters together in the system a — 
{Hq, Hk,Mi, po}. Two systems a and a are called equiva- 
lent if they are indistinguishable by all input-output 
experiments. Therefore by definition, we can estimate 
the real system a up to equivalence. Let us call the es- 
timated system a, which consists of estimated compo- 
nents a = {Ho,Hk,Mi,po}. We assume that the esti- 
mated system has been chosen to be of minimal dimen- 
sion, which implies that this system is also controllable. 
The goal is now to find a mathematical description of 
how different a can be from the real system a. 

Equivalence and similarity: — We first have to find a 
mathematical description of the equivalence. For some 
fixed input, equivalence means that the real system and 
the estimated system have to agree on all observable out- 
puts for all times, i.e., 

tT{Mip{t)} = tT{Mip{t)}, (1) 



where p{t) is the state evolving from the initial state po 
with the Hamiltonians Hq and Hk- This is not very 
useful mathematically, because it still involves solving 
the Schrodinger equation. There is an algebraic descrip- 
tion of this property that is much easier. Let us denote 
Ca = ■ ■ ■ , whcrc CK is a multi-index of length L 
with entries aj S 0, . . . , iV^. Further, we include the case 
i = as the identity superoperator and introduce sim- 
ilar notation Ca for the estimated system. Equivalence 
can then be formulated as 

tr{MiCaPo} = tT{MiCaPo} (2) 

for any sequence of the indices a. This can be thought 
of as an "infinitesimal version" of ([1]), and a simple proof 
of this statement is found in [3]. 

We call systems similar if and only if there is a simi- 
larity transformation between them 

Ck=TCkT-\ Me^MiT-\ Po^Tpo, (3) 

(fc ~ 0, . . . , Ni) where Aii and Aii represent the actions 
of Ml and Mg in the Liouville space. It is obvious that 
similarity implies equivalence. Similarity is much easier 
to handle than equivalence, because of its simple mathe- 
matical structure. 

Translation to quantum case: — In bilinear system 
theory it was shown that if a is controllable, then 
equivalence implies similarity. This is proven by explicit 
construction of the similarity transformation between a 
and a. Because there are some subtle differences in the 
quantum case, we briefiy repeat these arguments. 

Assume a and a are equivalent and pick an arbitrary 
state p. We show that due to the controllability of system 
(T, the state p can be expressed as 

p^^XaCaPo- (4) 

Firstly, because ipo E u(dimH) (the algebra of skew- 
hermitian matrices) the set R = {iA \ A — J2a ^aCaPa} 
is a subset of u(dim?^). Because = 
—i[[Hk, Hj], •] and we have controllability, the linear 
combinations of Ca include Cj^ = —i[H, • ] for any her- 
mitian H. This means that [iH, iA] e i?, so i? is an ideal. 
Because it is not equal to the identity and not su(dim'H), 
we must have R = u(dimH). Therefore, we can express 
any hermitian operator as XaCaPa, and in particular 
any state p, as in (|4]). 
We then define T by 

Tp^'^ XaCaPO- 
a 

There are many possible representations of p. In order 
to see that T is well-defined as a mapping, we need to 
verify that any two equal representations XaCaPo = 
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J2a ^a-^aPo imply ^a^aPo = ^a^aPo- By lin- 
earity, it is enough to show that 

>^aCaPO = =^ ^ XaCaPo = 0. (5) 

Suppose that '^a^aPo — 0. Then, we have for any f3 

tr< MlCf3 ^ XaCaPO \ = ^^^a tr{MeCf3aPo} = 0, 
^ Q! a. 

where Cpa = ^/3^a- Now, due to the input-output 
equivalence ^ between the two systems a and a, i.e., 
tT{MiCi3aPo} = trlMiCpaPo} for any /3a, we get 

tr<^ MiC/3 ^ XaCaPo \ = '^^a tT{Me£/3aPo} = 0. 

Since this holds for any (3 and the system is controllable, 
we conclude J2a ^a^aPa = 0, which completes the proof 
of ([5]). The mapping is onto due to the controllability of 
the system, and is shown to be one-to-one by reversing 
the argument which proved that it is well defined. Finally 
using controllability and the property Tpo = po it is easy 
to see that T has to fulfill ^ . 

Unitarity: — Since we restrict ourselves to unitary 
dynamics, it is possible to prove that the above similarity 
5( • ) = T • T^^ is actually inducing a unitary transfor- 
mation on Hamiltonians. First, we note that both the 
real and the estimated Liouvillians have the commuta- 
tor structure Ck = —i[Hk,»] and Ck = —i[Hk,'], be- 
cause we restrict ourselves to unitary dynamics. These 
Liouvillians form a subspace U of all possible Liouvil- 
lians. We first show that controllability implies that this 
subspace is mapped into itself by the similarity transfor- 
mation S. Indeed, a simple expansion of commutators 
combined with S being a similarity transformation shows 
that 

= -i[[Hk,H,],*]eU. (6) 

By linearity, any element H of the generated algebra has 
the property that S{—i[H,»]) G U. Because the system 
is controllable, this algebra is just the set of all trace- 
less hermitian matrices, and therefore S{U) = U. Since 
there is an isomorphism between Cfj = — and i?, 
we can represent the action of 5 on Z// by a correspond- 
ing action on su(dim'H). By linearity, this must be a 
linear and invertible map S. Indeed, from ([BJ it follows 
that S{[H,H']) = [S(H),S{H')]: S* is a Lie automor- 
phism. A theorem in [10[ states that all automorphisms 
on gl{n) (the general matrix algebra) are of the form 
S{X) = AXA-^ or S(X) = -AX^A'^. Our automor- 
phism is instead on the sub-algebra su(dim "H). By choos- 
ing a hermitian basis of gl{n) we can extend it uniquely to 



one of gl{n) and apply the theorem. The additional her- 
mitian structure demands furthermore that A~^ = A'^ . 
Thus S{H) = UHW or S{H) = -UH'^W. The latter is 
excluded because it would not preserve the trace of quan- 
tum states. Hence, under the premise of controllability, 
two systems are indistinguishable if and only if they are 
related through a unitary transformation 

Hk - UHkU\ Me - UMiU\ po - UpoU^. 

Usage of a priori knowledge: — In practice, it is rea- 
sonable to assume that some elements of the black box 
are known. Each known element shrinks the set of 
possible unitary transformations, because, for example, 
Hk = UHkW = Hk implies [C/, Hk] = 0. 

As an example, we consider two qubits coupled by an 
unknown Hamiltonian. We estimate them by performing 
arbitrary operations Hi = Xi® I2 and H2 = Yi® I2 on 
the first qubit and by measuring a) Zi(E)t2, b) Zi(S}Z2, or 
c) |*~)i2(*~|, where Xi, Yi, and Zi are the Pauli opera- 
tors of qubit i = 1,2, and |*~)i2 = (|01)i2 - |10)i2)/v/2 
is the singlet state. Assuming that the system is control- 
lable, we can apply the above results. 

First, the conditions [U, Hk] = reduce the unitary 
transformation U to li (g) J72j where U2 is a unitary op- 
erator acting on the second qubit, which may be param- 
eterized as U2 = e~5^"' '^2 with a unit vector n. We then 
impose another condition [U, M] = 0: 

a) In the first case with M = Zi ® I2, this condition is 
already satisfied and the unitary transformation U is not 
reduced any further, U = liigje"^^" '^^, where remain 
three parameters. 

b) In the second case with M = Zi(E) Z2, the condition 
reduces C/ to li ig) e~^^^^ with a single parameter. 

c) Finally, in the third case with M = we 
have 

[[/, M] = z sin ^ (n. |*+) 1 - n,\<p-){-i>- \ 

+h.c., 

which is vanishing only when sin(6'/2) = 0, i.e., U = 
1 up to an irrelevant phase. This shows that the Bell 
measurement is more efficient to estimate the system. 

Infection criterion for arbitrary systems: — Let us 
consider another more general example, a generic Hamil- 
tonian of a d-dimensional Hilbert space in the form 

Ho^ ^ c„rn\n){m\, (7) 

(Ti,m)eE 

where the orthonormal basis |7i) may be thought of as 
"local," and E are the edges of the graph G = {\n),E), 
that describes the non-zero off-diagonal {n 7^ m) cou- 
plings Cnm- We assume that a set of nodes C can be 
controlled {Hk = \k){k\, k € C), and that (at least) 
one particular node of C can be measured (Mi = |1)(1|, 
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1 € C). The crucial assumption about the set C is that 
it is "infecting" G ■ This property is defined by the fol- 
lowing propagation rules: 1) C is "infected"; 2) infected 
nodes remain infected; and 3) the infection propagates 
from an infected node to a "healthy" neighbor iff it is its 
only healthy neighbor. For an arbitrary Hamiltonian we 
can always find an infecting set; how many nodes it con- 
tains depends on how sparse the Hamiltonian is in the 
particular basis of consideration. In practice there are 
physical choices of the basis corresponding to local oper- 
ations, and many Hamiltonians are infected by acting on 
a vanishing fraction of nodes only. 

Based on the assumption that C is infecting one finds 
that the system is controllable, so our theorem can be 
applied. Firstly, there is a fc e C that has a unique 
neighbor I outside C. For that k we have [iHk,iHo] = 
-J2m&n{k)i'^kTn\k){m\ - c„ifc|m)(fc|), where n{k) is the 
neighborhood of k. Commuting it with iHk again yields 

\k){m\ + \m){k\). 

mGn{k) 

(8) 

For m e n{k) f] C\ on the other hand, we can single out 
terms by 

[[iHrn,iHQ],iHk] = ~i{ckrn\k){m\ +C„fc|TO)(fc|). 

By adding these to ^ for all m G n(fe) n C, only a single 
term i{cki\k){£\ + Qfc|^)(fc|) is left. Commuting this with 
iTJfc again gives Cke\k){£\ — C(k\i){k\. Finally, commuting 
the latter two and subtracting the term proportional to 
iHk we are left with i\£){£\. By induction, we can obtain 
|n)(n|, Vn. This implies full controllability 

If we assume that Hk and Mi are known, we need to 
look at the unitaries that commute with these operators. 
There will be many. However, we will assume here the 
knowledge that the Hamiltonian Hq has the form given 
in ([7]). Hence, we are talking about an indirect coupling 
strength estimation , where the topology E is known 
while the parameters are unknown. Let us see what this 
knowledge implies. Firstly, we have to have [Hk,U] — 
= —[Hk,U^] {k G C). Since Hk are projectors that 
implies that \k) must be an eigenstate of U and [/^ for 
all k d C. The estimated Hamiltonian Hq = UHqU^ 
has to be of the form Hq = J2{n,m)eE Cnm\n){m\, where 
Cnm are unequal to zero and could in principle differ from 
Cnm- The edges E must be the same for both Hq and Ho 
because we assume knowledge of the topology. Because 
C is an infecting set, there is one k G C that has a 
unique neighbor £ outside of C. The corresponding term 
in the Hamiltonian Hq is Cki\k){£\ + cik\£){k\. Because 
|fc) is an eigenstate of U this transforms under U • W 
into Ckie"'"'\k){£\W + cike-''^''U\£){k\. Because the edges 
E are the same for Hq and Hq there is a corresponding 
term Cki\k){£\ + cik\£){k\ in Hq. Furthermore, since \k) 



is an eigenstate of W no other node \n) can be brought 
to \k), i.e., {k\U\n) = 0. Given that £ is the only node 
outside C coupled to k we conclude 

Ckie'''-\k){£\U^ + Qfce-*-C/K)(fc| = Ck,\k) {£\ + cik\£) {kl 

which implies that \£) is an eigenstate of U. Finally, by 
induction we get that U must be a diagonal matrix in the 
"local" basis |n). Thus, up to the local phases of the basis 
vectors the Hamiltonian Hq is uniquely estimated. What 
is remarkable here is that we do not have to assume the 
knowledge of the phases of Cnm and it suffices to measure 
a single node. This generalizes the previous results 
substantially. 

Conclusion: — We have shown that controlable closed 
quantum systems can be estimated, in principle, up to 
unitary conjugation. This provides an easy-to-check cri- 
terion for experimental setups, telling us which extra con- 
trols, measurements or a priori knowledge are needed to 
achieve highly accurate quantum system identification re- 
quired for quantum computation. We have applied this 
criterion to a simple setup to show how the structure of 
the measurement observables can be important for esti- 
mation efficieny, and we have constructed a method to 
obtain a fully controllable and estimable system from an 
arbitray Hamiltonian. 
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